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SUMMARY 


This  report  describes  the  re-entry  physics  research 
activities  at  the  Heliodyne  Corporation  under  Contract  DA-04-49b- 
AMC-4b8(Z}  during  the  time  interval  of  12  June  to  31  December 
19b4.  Investigations  were  conducted  in  the  following  areas: 
mixing  and  chemical  reactions  in  turbulent  wakes;  stability  of 
two-dimensional  and  axi-syn  metric  laminar  wakes;  the  influence 
of  non-adiabatic  collisions  in  the  vibrational  and  dissociativ'e 
relaxation  of  a  high  temperature  gas;  analyses  of  optical  scintill¬ 
ation  data  from  re-entry  field  tests;  and  shock  tube  investigations 
of  the  spectroscopic  and  electromagnetic  properties  of  ablation 
product/air  mixutres. 

The  individual  research  topics  are  described  in  the  follow¬ 
ing  pages  in  the  form  of  abstracts  from  completed  reports,  or 
progress  reports  on  continuing  activities. 
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mixing  and  fluctuations  in  . 


A  TURDULENT  wake’ 


by 


Andrew  P.  Freudian 


he  structure  of  turbulent  wakes  of  hypersonic  b  d' 
argued  to  be  lat--,  i  bodies  is 

O  be  largely  controlled  by  the  combined  actio,  r  , 

large  scale  action  of  relatively 

rirrrdTT'  — 

--lent  diatortmna  oI  Ilmd  eleme:: 

-thematical  model  representing  I^rbulent TalrmlTT^' 

a  mixing  lag  is  discussed  ,4  ‘  ®  terms  of 

K  aiscussed  and  qualitatively  iustifi^d  • 
above  described  •  •  ^  justified  m  terms  of  the 

M  ’  ^^^^e'vhat  generalized 

Numerical  predictions  for  the  ma.n'i  ^  ieneraliaed. 

icr  me  magnitude  and  variation  ^ 

density  fluctuations  in  a  relative-v  lo  . 

.  1  •  eiative.y  low  speed  wake  are  obt;,:raed  v 

-Giving  a  simplified  set  of  obtained  by 

P  ed  set  of  equations,  in  which  the  lay  ic 

dy  an  effective  mixing  boundary,  and  are  in  -P-^-‘ed 

-‘H  the  general  magnitude  and  trend  o 

although  firm  experimental  mea  ^Perimentally, 

h-ci  i.iiiencai  measurements  i  i  • 

menis  are  lacking  to-date. 


— 

Abstract  o*  Heliodyne 


Research  Report  8, 


dated  December  I964 
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ANALYSIS  OF  PARTIAL  WAKPw  MIXING  AMD  APPLICATION 
TO  REACTING  NON  EQUILIBRIUM  WAKES" 


by 

A,  Proudian  and  A.  Widawsky 

The  description  of  the  mean  properties  of  turbulent  hyper¬ 
sonic  non-equiUbriurr>  wakes  in  terms  of  partial  mixing  of  ingested 
fluid  by  the  growing  turbulent  core  is  analyzed.  A  formulation  of 
the  conservation  equations  for  wake  mixing  is  developed  which 
avoids  internal  inconsistencies  in  the  specification  of  each  fluid 
portion  present  in  some  existing  treatments.  The  problem  of 
representing  the  important  effects  of  statistical  fluctuations  in  mix¬ 
ing  in  the  conservation  equations  for  the  mean  wake  properties  is 
briefly  discussed.  The  conservation  equations  in  the  case  of  n«^gli- 
gible  dynamic  interaction  between  the  mixed  and  unmixed  fluids  and 
in  the  limit  of  si  dden  mixing  are  integrated  to  determine  mean  wake 
properties,  particularly  electron  densities  and  their  fluctuations, 
for  a  blunt  body  case  of  importance.  The  sudden  mixing  limit  is 
argued  to  be  a  reasonable  first  approximation  to  turbulent  wake  mix¬ 
ing.  The  application  of  the  results  to  the  prediction  of  radar  back- 
scatter  from  turbulent  underdense  wakes  is  discussed. 


Abstract  of  Hcliodyne  Research  Note  16, 


dated  December  1964 
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ANALYSIS  OF  SCHLIEREN  MEASUREMENTS^^ 

by 


Andrew  Proudian 


7-cu.e;.:::!''xbr‘t:L‘T.;.. 

bo  dominant  in  the  oL'cn  “d  ‘"'Pd^tant  and  may  even 

cast  some  doubt  on  the  vaUdUrof  th  \  appears  to 

Clav  e.  »l  „f  ,  V  ‘  ‘"‘arpretations  placed  by 

Clay,  etal,  of  Lincoln  Laboratory  on  their  ob.  erved  results  A 

"  or'd^rt  -Reaction  r  ;<,uted 

in  order  to  arrive  at  a  definite  conclusion.  requ.rcc 


^Abstract  of  Heliodyne  Res 


earcli  Memorandum  8, 


dated  October  1964 
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STABILITY  OF  DUST  LADEN  TWO-DIMENSIONAL 
LAMINAR  WAKES'" 

by 

All  Hasan  Nayfeh 

The  stability  of  two-dimensional  compressible  dusty  wakes 
is  investigated.  An  equation  is  derived  for  the  amplitude  of  the 
pressure  disturbance.  The  solution  of  this  equation  is  sought  by 
perturbing  the  dust-free  solution  for  small  mass  fraction  of  dust  to 
gas.  The  straightforward  perturbation  expansion  is  not  uniformly 
valid  far  away  from  the  wake  axis.  A  second -order  uniformly  x^lid 
solution  is  obtained.  The  results  show  that  the  dust  can  be  stabi¬ 
lizing  as  well  as  destabilizing,  depending  on  the  size  of  the  particles 
and  their  distribution.  For  a  uniform  distribution  of  dust  and  for 
%''ery  small  particle  size,  the  dust  has  a  destabilizing  effect.  As  the 
dust  size  increases,  the  destabilizing  effect  increases  to  a  maxi- 
m.um,  and  then  decreases  to  zero.  As  the  particle  size  increases 
further,  the  dust  will  have  a  stabilizing  effect.  The  quaiititati ve 
stabilizing  or  destabilizing  effect  of  the  dust  depends  on  the  velocity 
defect  and  temperature  cxces=^  at  the  wake  axis. 


Abstract  of  Heliodyne  Research  Report  6,  dated  October  1964 


b 
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STABILITY  OF  DUST  LADEN  AXI-SYMMETRIC 
LAMINAR  WAKES* 

by 

Ali  Hasan  Nayfeh 


in  order  to  investigate  the  stability  of  dusty  axi-symmetric 
laminar  wakes,  a  single  equation  for  the  amplitude  of  the  pressure 
turbance  is  derived.  The  stability  problem  is  converted  into 
finding  the  eigenvalues  (wave  number,  complex  wave  velocity)  and 
the  eigenfunctions  of  a  complicated  eigenvalue  problem.  Thus,  the 
aolution  is  sought  by  perturbing  the  dust-free  solution.  The  straight- 
forward  perturbation  expansion  breaks  down  far  away  from  the  wake 
axi-s.  A  uniformly  valid  second  order  expansion  is  obtained,  and  the 
characteristic  equation  for  the  eigenvalues  is  derived.  Thus,  the 
aiialytical  analysis  is  completed. 

For  a  Gaussian  mean  flow  model,  the  eigenvalues  and  ampli  - 

Xic3.tiori  rector  will  *-  t  3 

oe  ootained.-  The  effect  oi  dust  on  the  position  of 

the  transition  point  behind  a  body  will  be  investigated. 


if 

Summary  of  continuing  work 
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HIGH  '^MPERATURE  MOLECULAR  COLLISIONS 

by 

A.  Proudian  and  I.  Tannenbaum 

Calculations  of  the  rates  of  dissociation  and  vibrational 
transitions  in  gases  at  several  thousand  degrees  Kelvin  have  been 
undertaken  in  the  preceding  period.  The  objective  of  the  calcula¬ 
tions  has  been  primarily  to  include  the  effects  of  the  higher  energy 
molecular  collisions  which  become  important  at  high  temperatures 
in  determining  the  population  of  the  various  molecular  vibrational 
levels,  and  indirectly  affect  the  state  of  a  high  temperature  non¬ 
equilibrium  gas.  Those  collisions  particularly  affect  the  population 
of  the  higher  vibrational  levels,  which  dominate  the  dissociation 
process.  Previous  calculations  of  coupled  dissociation  and  vibra  • 
tional  relaxation  have  used  models  of  molecuL!,r  collisions  which  are 
largely  adiabatic  in  nature,  and  which  do  not  strongly  allow  multiple 
quantum  transitions,  in  contrast  to  the  present  model  which  assumes 
♦Lat  impvilsive  collisions  dominate  the  relaxation  process,  particu¬ 
lar,  ly  insofar  as  the  upper  vibrational  levels  are  concerned. 

The  study  of  the  coxipled  dissociation-vibrational  relaxation 
problem  requires  the  determination  of  transition  probabilities  in 
impulsive  collisions,  together  with  the  solution  of  a  set  of  so-called 
Master  Equations  describing  the  relaxation  process.  The  present 
approach  to  the  computation  of  the  transition  probabilities  has  been 
to  assume  that  the  inter -molecular  forces  are  large  compared  to  the 
molecular  binding  forces  during  the  molecular  collisions.  It  has 


Summary  of  continuing  work 
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also  been  assumed,  (consistently,  )  that  interaction  occurs  between 
on.y  one  nucleus  of  each  member  of  a  colliding  pair  of  molecules, 
and  that  the  interaction  potential  is  independent  of  the  position  of  the 
non  colliding"  nuclei  of  the  colliding  molecules.  Thus  it  has  been 
umed  that  multiple  collisions  are  not  important.  The  possibility 
ol  collision  of  either  member  of  a  molecular  pair  with  either  mem¬ 
ber  of  the  collision  partner  is  then  accounted  for  by  simply  adding 
the  inditi„ual  cross  sections  for  each  possible  collision.  The 
probability  of  transition  from  a  given  state  i  to  a  state  f  of  one  of 
the  colliding  molecules  can  be  equated  within  the  validity  of  the 
above  stated  assumptions,  to  the  probability  of  transition  between 
free  state  wave  packets  with  momentum  distributions  corresponding 
to  the  bound  states  i  and  f  integrated  over  the  range  of  momenta 
represented  by  the  bound  state.  Collisions  in  which  both  molecules 
undergo  an  internal  change  are  assumed  to  have  a  imail  cross  sec¬ 
tion  and  are  neglected,  at  least  for  the  present.  At  present,  rota- 
tionai  transitions  are  also  ignored,  and  the  bound  state  wave  func¬ 
tions  to  be  used  are  based  on  a  Morse  intranuclear  potential.  On'y 
qua.si-hard  sphere  collisions  have  been  considered  to  date,  as  the 
free  state  scattering  amplitudes  are  then  known  analytically.  The 
nteraction  is  denoted  as  quasi-hard  sphere  rather  than  simple 
hard  sphere  because  the  radius  of  the  hard  sphere  is  taken  to  be  the 
olasstcai  turning  radius  for  motion  in  a  Morse  or  Leonard  Jones 
interaction  potential.  Thus  the  hard  -phere  radius  is  made  to 
depend  on  the  incident  energy.  This  should  provide  a  substantial 
imp-  ,  -ement  over  a  simple,  fixed  radius  hard  sphere  model.  A 
more  realistic  potential  with  adjustable  parameters  will  be  dev.  oped 
for  use  ,n  applications  to  specific  molecules  of  interest.  An  optimum 
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numerical  method  using  the  hard  sphere  amplitudes  for  computing 
the  transition  probabilities  is  currently  being  sought. 

The  current  formulation  of  the  relaxation  equations  has 
been  made  to  correspond  to  the  relaxation  of  a  polyatomic  gas  behind 
a  normal  shock.  At  present  Uie  formulation  applies  to  a  system  of 
molecules  with  only  one  internal  set  of  energy  levels,  e.  g.  ,  the 
vibrational  levels  of  a  diatomic  molecule.  However,  the  formula¬ 
tion  ca*.  be  readily  generalized  to  include  more  tha  one  vibrational 
mode,  or  rotational  modes. 

The  formulation  applies  equally  to  molecules  in  an  infinite 
heat  bath  at  a  fixed  temperature  T  Trid  to  an  isolated  system  of 
molecules  for  which  T  will  drop  as  the  gas  relaxes. 

The  key  simplifying  assumption  of  previous  calculations  of 
coupled  disGoci"  ion  and  vibrational  relaxation  haG  been  that  a  known 
form  for  distribution  of  the  molecules  among  their  vibrational  states, 
generally  a  Boltzmann  distribution  corresponding  to  some  tempera¬ 
ture  T  ,  was  assumed,  and  therefore  a  single  parameter,  such  as 

V 

T  was  solved  for  in  order  to  obtain  the  vibrational  populations. 

This  simplification  is  not  consistent  with  the  present  model  for 
transition  probabilities  and,  in  fact,  the  study  of  the  actual  distri¬ 
bution  is  one  objective  of  this  study.  Thus,  the  population  of  each 
vibrational  level  must  be  considered  individually,  leading  to  a  set 
of  up  to  fifty  or  more  coupled  differential  equations  to  describe  the 
system.  Since  these  cannot  be  integrated  analytically,  care  must 
be  exercised  n  selecting  a  numerical  integration  scheme  which  is 
not  prohibitively  time  consuming.  At  present,  the  numerical  inte¬ 
gration  of  the  relaxation  equations  is  being  studied  for  the  simplest 


12-64 


9 


c«se,  which  corcesponda  ic  coupled  ■.bratioo-  J.asocia.ion  relaxa- 
Uon  Of  a  highly  diluted  diatotoic  gas  an  infinite  ntonatotnic  (pec- 
ect)  heat  bath.  The  initial  computations  will  be  performed  for  an 
oxygen  argon  mixture,  for  which  experimental  results  are  available. 

The  dissociation  rates  to  be  used  in  solving  the  relaxation 
equations  have  not  yet  been  determined.  While  the  impulsive  colli- 
Sion  model  can  be  applied  to  predict  dissociation  probabilities,  the 
criteria  for  its  validity  are  not  well  satisfied  for  dissociation  from 
most  of  the  vibrational  levels  except  for  collision  energies  of  per¬ 
haps  ten  volts  or  so.  When  numerical  results  from  the  collision 
model  for  various  temperatures  are  available,  it  will  be  possible  to 
choose  between  using  those  results  and  relying  on  experiments 
equilibrium  rates  together  with  a  preferential  dissociation  model 
such  as  that  used  by  Treanor  and  Marrone.  which  favors  dissocia¬ 
tion  from  higher  vibrational  levels,  and  which  is  consistent  with 
theoretical  numerical  results  such  as  those  of  Keck. 

It  is  hoped  to  compare  the  resuUs 

t'  rebuacs  with  the  experimental 

studies  by  K.  Wray  of  relaxation  of  dilute  oxygen  in  argon  behind 
shocks,  and  with  the  theoretical  predictions  based  on  the  Landau- 
teller  transition  probabilities,  if  the  numerical  calculations  for 
the  simple  cases  indicate,  as  anticipated,  a  significant  change  in 
dissociation  rates,  the  formulation  will  be  ganeraliaed  to  flow  situa¬ 
tions  of  interest  behind  hypersonic  shocks  and  in  hypersonic  bound¬ 
ary  layers,  and  applmd  to  ablation  contaminated  flows  and,  in  par- 

ticular,  to  the  study  of  the  dissociAfirxn  e  ■ 

y  oi  me  aissociation  r  important  radiating  states 

of  excited  molecules. 
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OPTICAL  SCINTILLATION  ANALYSES'!' 


by 


Saul  Feldman 


Analyses  relating  the  time  dependence  of  passively  emitted 
radiation  to  significant  body  parameters  of  a  re-entering  vehicle 
have  recently  been  initiated,  anci  no  definite  conclusions  hove  yet 
been  deduced. 


-Summary  of  continuing  work 
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A  SUMMARY  OF  CALCU  LATION3  OF  CHEMICAL  AND  FLUID 
DYNAMIC  PARAMETERS  IN  THE  SAP  AG  FACILIIY^:^ 

by 

William  J.  Hooker 

The  equations  describing  the  ablation  lifetime  and  trajectory 
material  particles  in  a  shock  tube  flow  (developed  in  previous 
reports  on  the  SAP  AG  facility)  have  been  rc -examined,  and 
programmed  for  computation  on  an  IBM  7094.  Computations 
have  been  made  of  particle  burnup  times  for  various  materials 

e^sL^k"'  b  =  compared  with  calculat- 

het  °  tarries.  In  the  limit  of  small  relative  velocity 

between  the  particles  and  the  free  stream,  simple  analytic  sol- 
u  ions  exist  for  particle  radius  and  trajectory  which  are  applicable 
to  many  physically  useful  shock  tube  conditions. 


^Abstract  of  Heliodyne  Research  Note  17,  dated 


January  1965 
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EXPERIMENTAL  MEASUREMENTS  IN  THE  SAPAG  FACILITY*-’* 

by 

W.  J,  Hooker,  A.  L.  Morsell  and  1.  R.  Tannenbaum 

The  final  version  of  the  SAPAG  facility  has  been  in  operation 
since  October  1964,  The  basic  shock  tube  has  an  all-stainless 
steel  driven  section  of  3"  diameter  x  30'  lengtn,  and  is  terminated 
with  a  dump  tank.  The  tube  currently  is  cold  gas  driven,  and 
shock  waves  with  velocities  up  to  4.  Z  mm^sec  have  been  generated 
in  air  using  a  hydrogen  driver  at  100  atmospheres  pressure. 

The  problem  associated  with  producing  powders  in  the 
submicren  size  range,  storing  them  for  extended  periods  of  time, 
and  injecting  them  into  the  showk  tube  in  a  uniform  and  reproducible 
manner  have  basically  been  overcome.  At  the  present  time,  powd¬ 
ers  are  being  produce-i  that  have  a  particle  size  distribution  such 
that  approximately  94%  of  the  par*^icles  are  smaller  than  0.  8  micron 
diameter  and  50%  smaller  than  0 . 1  micron,  as  determined  by- 
electron  microscopy.  Storing  of  the  finely  ground  powders  for 
periods  of  time  up  to  six  weeks  produces  no  significant  change  in 
the  particle  size  distribution. 

tVhereas  previous  tests  of  powder  injection  uniformity  in  a 
gless  shock  tube  produced  substantial  gradients  of  powder  mass 
loading  in  the  axial  direction,  measurements  m  the  present  all 
metal  facility  sliow  less  than  a  10%  variation  from  ore  end  of  the 
shock  tube  to  the  other,  and  the  reproducibility  is  better  than  10%. 

Measurements  have  been  made  of  the  rate  of  adsorption  of 
water  vapor  into  the  surfaces  of  the  finely  ground  pow-ders  used  in 
experiments  to  date.  These  tests  have  shown  that  for  particles 
prepared  under  drybox  conoitions,  there  is  a  rapid  adsorption  in 
the  first  few  seconds  of  exposure  to  a  normal  room  atmosphere, 
but  that  the  final  equilibrium  weight  of  any  given  sample  corre¬ 
sponds  to  the  addition  of  approximately  one  molecular  layer  of 
'vater  vapor  to  each  particle.  For  the  materials  and  sizes  tested 
to  date,  this  corresponds  to  a  w'eight  increase  ot  approximately  1%. 

*l*Sv^mmary  of  continuing  work. 
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Extensive  investigations  have  been  made  of  shock  wave 
attenuation  along  the  length  of  the  tube,  and  testing  time  at  the 
observation  stations.  For  initial  pressures  of  1  to  ^00  mm  Hg 
and  for  shock  velocities  from  1.  8  to  4.  2  mm^sec,  the  attenuation 
IS  less  than  b,,  over  the  useable  length  of  the  tube  (aoproximately 
e  last  half)  The  measured  testing  times  fall  approximately  half 

way  e  ween  t  e  calculated  values  based  upon  laminar  and  turbulent 
boundary  layer  flows. 

shock  tube  has  been  instrumented  with  a  Jarrell-Ash 
/6.  3  plate  spectrograph,  and  an  EG&G  xenon  flash  lamp  for  ab¬ 
sorption  measurements.  The  flash  lamp  has  a  nominal  1/e  dur- 
a  ion  of  l^sec  and  is  triggered  to  flash  at  a  selected  time  behind 
the  incident  shock  wave.  Emission  and  absorption  spectra  have 
been  obtained  for  shock  heated  air  and  teflon,  and  air  and  epoxy 
resin  mixtures,  using  a  10  micron  spectrograph  slit  v/idth  (with 
a  basic  grating  resolution  of  20A/mm).  Photomultipliers  and 
intrared  detectors  with  appropriate  wavelength  interval  filters  are 
currently  being  selected  to  follow  reactions  pertinent  to  each 
material  being  tested. 

The  basic  circuitry  and  configuration  of  the  coil  used  for 
electron  density  measurements  in  the  shock  tube  have  undergone 
extensive  revision.  Various  coil  geometries  and  circuit  elements 
have  been  tested  from  the  point  of  view  of  basic  circuit  stability 
as  well  as  the  magnitude  of  the  frequency  shift  under  actual  shock 
tuoe  conditions.  At  the  present  time,  no  single  configuration  works 
ccep  aoly  over  the  entire  range  of  electron  densities  encountered 
(based  upon  pure  air  equilibrium  chemistry),  and  it  may  be  neces¬ 
sary  to  use  selected  designs  for  several  overlapping  ranges  of 
operation. 


A  series  of  test  firings  into  pure  air  have  been  run  in  order 
to  calibrate  the  electron  density  measuring  coil  under  the  proper 
conditions  of  electron  density  and  collision  frequency.  From  the 
tests,  a  curve  of  coil  frequency  change  versus  calculated  electron 
ensity  is  obtained.  Since  the  data  interpretation  in  this  calibration 
technique  requires  an  equilibrium  electron  density  to  be  establish- 
e  during  the  available  test  time,  a  nonequilibrium  normal  shock 
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code  obtained  from  Cornell  Aeronautical  Laboratory,  Inc.  (CAL 
Report  No.  QM-1626 - A-12(II)  "Inviscid,  Nonequilibrium  Flow 
Behind  Bow  and  Normal  Shock  Waves,  "  by  Leonard  J.  Garr  and 
Paul  V.  Marron*  )  has  been  used  to  compute  the  equilibration 
distance  for  the  test  cases  of  interest.  It  is  anticipated  that  this 
computer  code,  along  with  the  normal  shock  code  obtained  from 
the  Los  Alamos  Scientific  Laboratory  (Report  No.  LA-2980,  "HUG, 
A  Fortran-FAP  Code  for  Computing  Normal  Shock  and  Detonation 
Wave  Parameters  in  Gases,  "  by  Paul  F.  Bird,  Russell  E.  Duff 
and  Garry  L.  Schott),  will  be  extensively  used  in  this  measurement 
program. 


